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Abstract: The relatively clean theoretical probes of the Standard Model (SM), and the 
various theories beyond the SM, provided by radiative, semi-leptonic and (purely) leptonic 
decays of B-mesons have become increasingly important. Due to the large number of 
possible distributions in the semi-leptonic decays based on the quark level transition b — > 
sl'^£~ (not just the branching ratio), these transitions have become very useful. A study 
of the Forward-Backward asymmetries for the inclusive decay {B — > Xsl~^l~) is carried 
out in this paper. This study shall be performed in the SM and a minimal supersymmetric 
extension of the SM, namely the mSUGRA model. 
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1. Introduction 

The recent flux of papers studying the various decay processes of B- mesons stands testimony 
to the great importance of this field in testing not only the parameters of the Standard 
Model (SM) but in constraining the parameters of its many possible extensions. In this 
vain processes of particular interest are the so called "rare" decays, these being transitions 
through loop levels only [1]. Of particular interest are the Flavour Changing Neutral 
Currents (FCNC) b s{d). These rare decays are extremely sensitive to the masses 
and couplings of the virtual particles involved in this transition, and therefore provide the 
greatest hope of determining many of the SM's (and its possible extensions) underlying 
structures. We therefore require a decay process which is both theoretically clean and 
experimentally realizable in the near future, such as B ^ X^-f^ and B XstT . Of 
these two processes this paper shall focus on the inclusive decay mode, B Xgi'^i", for 
the reason that this decay mode is a far richer alternative to the radiative mode as here 
many distributions involving the final state lepton pair can be measured. In the inclusive 
decay channel {B — > Xsi'^i~) observables such as Forward-Backward (FB) asymmetries 
and lepton polarization asymmetries, where one or both the leptons are polarized, have 
been extensively studied both within the SM and its various extensions [2, 4, 5]. 

Recently, it has been pointed out by Bensalem et al. [6] that one can construct the 
polarized FB asymmetries also (FB asymmetries when one or both the leptons are polar- 
ized). These new observables could provide a sufficiently large number of observables for a 
very strict testing of the SM and any new physics. These polarized FB asymmetries were 
studied in an exclusive decay mode {B Ki~^i~) both within the SM and the Minimal 
Supersymmetric SM (MSSM) [7], motivated by the vanishing unpolarized FB asymmetry 

^this decay mode has already been observed in various B-factories and has proved to be extremely useful 
in placing stringent bounds on many new physics models 



- 1 - 



for B — > Ki^£~ within the SM [8,9]. Based on the construction of doubly polarized ob- 
servables by Bensalem et al. such observables were studied in the case of i? — > K*i~^i~ [10] 
and in the inclusive mode [5]. In this work we will focus on the polarized FB asymmetries. 

Lately there has been an increased interest in how the predictions of various observ- 
ables will change in the two Higgs doublet model (2HDM) and Supersymmetric (SUSY) 
extensions of the SM. The primary reason for this is due to the increase in the number 
of particles in these theories, such as the Neutral Higgs Bosons (NHBs), and hence scalar 
and pseudo-scalar exchange operators are now included. Such operators give rise to many 
interesting new possibilities, like orders of enhancement in the branching ratio for purely 
leptonic decays {Bg^d £^l^) [11, 12] and non- vanishing values for unpolarized FB asym- 
metries in B ^ K('jr)l^£^ [7~9]. As the coupling of the NHBs to the leptons is proportional 
to the lepton mass [2, 11, 12, 14], the effects of these new particles is greater when the final 
state leptons are either ij, or r. Therefore in this paper we will also consider the effect 
of the new scalar and pseudo-scalar operators which arise from a MSSM. The particular 
model to be employed shall be explained in more detail in section ^. 

As such, this paper shall therefore be organized as follows: In section |2| we shall review 
the effective Hamiltonian and the matrix element for the process concerned, as derived 
from the quark level transition b s£~^£~. Section |^ will then focus on the definitions and 
analytic results of the FB asymmetries. Finally, section |^ will contain our discussion and 
numerical analysis of the results. 



2. Effective Hamiltonian 

The effective Hamiltonian for the quark level transition b s£~^£~, as first studied by 
Grinstein et al. (and later to NLO by Buras and Miinz) [15]^, and represents the quark 
level process for the inclusive decay B Xs£~^£~ , under consideration (where we shall use 
descriptions as penned in references [2, 11, 14]); 

/ 10 10 \ 

Heff = -^ytbv;, [J2ao, + Y^^CQ^QA (2.1) 

where the Oj operators are the current-current (i=l, 2), penguin (i=3,. . . ,6), magnetic 
penguin (i=7,8) and semi-leptonic (i = 9,10) operators with the Q corresponding to the 
standard Wilson coefficients. The additional operators Qi (i = 1,. . . ,10) and their Wilson 
coefficients {Cq.) are related to the additional particles of the MSSM under consideration. 

Using the above effective Hamiltonian, and neglecting the s-quark mass, the matrix 
element for the quark level transition b s£~^£~ is; 

M = ^VtbV;A -2C'/f'^{sia,,q''PRb){h^£) + Cl^f{s^^^ 
v2vr q 



+Cio{sif,PLb){£ri^£) + CQ,{sPRb){££) + CQ,{sPRb){h 




^the complete NNLL calculation of B ^ Xsi'^i has recently been given in reference [16], however for 
our analysis we will use only the Leading Log results given in reference [15]. 
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In this expression q represents the momentum transferred to the lepton pair, given as 
q = p+ + where p+ and p- are the momentas of the l'^ and l~ particles respectively. 
Note also that the CKM factors have been denoted VtbV^g above, and that Pl,r = (lT75)/2. 

In this paper we will also include the long range effects which arise from the four 
quark operators l~ s\Oi\b) for i = 1, . . . , 6. This shall be done by absorbing them into 
the description of the Wilson coefficients; a prescription which has been used in many 
earlier works, see references [4,17]. These long distance effects typically arise from cc 
contributions, and are taken into account by adding a term to the C^^^ coefficient; 

Cr oc n V ^^^"(^ ^ ^-^^)^ltai (2.3) 

where we have used the same symbols and notation as in Kriiger and Sehgal [4]. In the 
above equation the phenomenological factor k is introduced to predict the correct branching 
ratio of Br{B — > J/^Xs Xsi~^i~). For our numerical analysis we will choose k to have 
a value of 2.3. 

Using the expression of matrix element in equation (^]^) we obtain the expression for 
the differential decay rate as; 



- T9^4^'^*^^*^' (1-s) yi-^A (2.4) 

with 

A = 4^^ (1 + ^) \C^/f\^ + (1 + 2^) (1 + ^) ICI^^^ 

+(1 - Srhj + 2s+ ^-^)\C^o\' + ^(-4mf + s)\CqA' + ^^I^q/ 

+12(1 + ^)Re{Cl"*C'j") + 6meRe{CloCQ,). (2.5) 

Using this expression of the invariant mass spectrum (including the scalar exchange effects) 
we will now analyze the various FB asymmetries in the next section. 

3. Polarized FB asymmetries 

As was mentioned in the introduction, the additional intricacies of the b — > s£~^£~ quark 
level transition requires many more experimental observables in processes such as i? — > 
Xsi^i^ in order to constrain all the parameters of a particular MSSM. With this in mind 
we shall now define the differential FB asymmetry as [17]; 

A(s) - C - r ^dz (3 1) 

Jq dsdz y_i dsdz 

If we do not sum over the spins of the outgoing leptons then the FB asymmetry will, in 
general, be a function of the spins of the final state leptons; and can be defined as 

As-',s+rs)= / ' dz- / ' dz. (3.2) 



dsdz J_i dsdz 
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As our eventual aim is to derive expressions to be searched for in experiments it is conve- 
nient to use the normahzed FB asymmetry. Normahzing the above (equation 3.2) definition 
by dividing by the total decay rate we get; 



dr/ds 



(3.3) 



We can split the FB asymmetries into their various polarization components and we will 
do this in analogy to the prescription given in Bensalem et al. [6]; 



Ai^S ,5 ) — A ~\~ A-iSj^ -\- A-iS^ A-ijSj^ Sj 



(3.4) 



where i,j = L,T,N are the longitudinal, transverse and normal components of the polar- 
ization^. Therefore, from equation ( ^.41 ) we can write the single polarized FB asymmetry 
as; 



At = Ais- 



I, s 



--j)+A{s- 
i) +A{s' 



I, s 



-A{s- ■- 
-4(s" = i, s 



I, s 



--j)-A{s- 
-i) - A{s~ 



-z, s 



-J) 

(3.5) 
(3.6) 



and the double polarized FB asymmetry as; 
Aij = ^(s^ = z, = j)—A{s~' =i,s^ = —j)—A{s~ 



-I, s 



j)+A{s- 



-hs'^ = -j)- 
(3.7) 

Using the expression of the matrix element as given in equation ( |2.2D we can get the 
expression of unpolarized FB asymmetry; 



A 



Am 



^2 



A 

X 



The analytical results of the polarized FB asymmetries are then; 



AT = ^ 



e//|2 



+ 1 



Ami 



\Cio\' + 6Re{C^"*q 



9 



-12^i?e(C7^^^*Cio) - 6m^i?e(C77^^^*CQ,) - %^Re{Cl^ ^* C^^) 
+3meRe{Cl"*CQ, 



A 



N 



2 (1 — s)rhi 

A 71 



Ami 



Re{Cl"*Cio) 



(3.9) 
(3.10) 

(3.11) 



^the convention used here is that where a repeated index appears it is summed over 
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Ami 
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A 
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-2Re{Cf^Cio) 
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(3.12) 
(3.13) 

(3.14) 



2 (1 — s)mi 



Ami 



A VI 

—Aln 
3 



im{ci^a 



eff^ 
9 ) 



e//|2 



+ |C. 



e//|2 



m, (2iie(C7^^^*CQj +i?e(C9^^^*C7Qj) i / 1 



A 

6 . ^e// 



{i2e(C9^^^*CQj + i2e(C7^^^*CQj} 

(3.15) 

(3.16) 



(3.17) 
(3.18) 

(3.19) 



Ami; 



-rfn[2^Im{ClQC, 



Im{C^/f*CQ,) + !I^Im{C*,oC^/f: 



Atl = —Alt 
Atn = Ant 
Att = —Ann 



(3.20) 

(3.21) 
(3.22) 
(3.23) 



where A is given in equation ( |2.5[ ). We will discuss the above obtained expressions of the 
various FB asymmetries and present our numerical analysis of the same in the next section. 



4. Numerical analysis and Conclusions 



As it is experimentally more useful we shall present our results in the form of average 
values^. Though, as is well known [27], we have modelled the inclusive decay B Xs£^£~ 
from a quark level transition b — > si'^i" using the Heavy Quark Expansion (HQE) in 
{1/mi)), where in this expansion we have used the first term of this expansion. However, 
this expansion breaks down when s — > 1. In fact for higher values of s one should consider 
the non-perturbative (1/m^) corrections as well [27]. As was also shown in [27] it is difficult 

*it may be possible that the averaged values of the observables is small but still there could be large 
predicted values of the same observables in certain dileptonic invariant mass regions 
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Figure 1: The variation of the branching ratio of _B ^ XgT^r^ with tanP in the mSUGRA model. 
Other model parameters are: TOi/2 = 450GeV, A — 200GeV and the sign of /i will be taken to be 
positive for all our numerical analysis 



to estimate these non-perturbative corrections for s > 0.8 . For our analysis we have placed 
a cutoff on the upper limit of s. The averaging procedure which we will be using is defined 
as: 



0.8 



(3.646+0.02)2/m| 



A^ds 

2 as 



0.8 



dV 



(4.1) 



(3.646+0.02)2/m2 ds 



ds 



that is, in the calculation of our averages we have taken the lower limit of integration to 
be above the first resonance^. Note that we have used the input parameters presented 
in appendix The results of table |l|, these being our SM predictions of the integrated 
observables, will of course be altered when considering a MSSM. 



Bv{B XsT+T-) 


A 


^1 




All 


Aln 


Alt 


Ann 


Ant 


1.8 X 10"'^ 


- 0.176 


0.455 


-0.029 


0.176 


-0.0084 


0.063 





-0.083 



Table 1: Our SM predictions of the integrated observables 

We have also performed the numerical analysis of the SUSY effects on various ob- 
servables, presented in the previous section, for the inclusive mode {B XsT^t^). If 
we assume that SUSY exists, then the MSSM is the simplest SUSY extension of the SM. 
However, the MSSM has a large number of parameters, which limits its practical usage for 
phenomenological studies. There are many models which reduce the vast MSSM parameter 
space to a manageable set of parameters. These include Dilaton, Moduli, mSUGRA (min- 
imal Supergravity), rSUGRA (relaxed SUGRA), CMSSM (constraint MSSM) etc. All of 
^by first resonance we mean the resonance after the threshold of the decay, which is s > 4m^. 
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Figure 2: Ann variation with the invariant mass of the dileptons, s, in the mSUGRA model. The 
other model parameters are: mo = 400GeV, mi/2 = 500GeV, A = 0. 

these models have their own characteristics. The basic feature of all these unified models is 
that they assume some unification of the parameters at a higher scale (usually at the GUT 
scale). In our numerical analysis we have used one of the more popular models, minimal 
Supergravity or mSUGRA^. 

In the mSUGRA model one assumes the universality of masses and coupling constants 
at the GUT scale. The mSUGRA framework has five independent parameters'' namely, 
mo (the unified mass of all the scalars), mi/2 (the unified mass of all the gauginos), A 
(the unified trilinear coupling constants), ian/3 (the ratio of vacuum expectation values of 
the two Higgs doublets). Using these unified parameters the renormalization group (RG) 
equations of all the parameters are evolved from the GUT scale to the electroweak scale and 
then at the electroweak scale one imposes the condition of electroweak symmetry breaking 
(EWSB). Imposition of EWSB fixes the magnitude of the parameter /x ®, but the sign of 
H still remains arbitrary. As such the sgn{fj,) ^ is taken to be another parameter, which 
makes for five parameters (out of these five, four are continuous parameters with the of 
sgn{^) being discrete, that is ±1). 

Furthermore, we have worked in the high tan[3 region of the mSUGRA parameter 
space. This was done as it is only in this region that the contributions of the Neutral Higgs 
Bosons (NHBs) become significant [2,10-12,14]. The parameter space of mSUGRA has 
been constrained by the experimental observation of i? — > Xg^f [19]. For our numerical 

^details of the model can be found in Nilles [22] 

^in Moduli and Dilaton scenarios one more additional condition is used to reduce the number of param- 
eters to four 

*this is the coupling parameter of the two Higgs doublets 

®the convention for sgn{fi) which we will be going to take is such that fi appears in chargino mass matrix 
with positive sign 
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Figure 3: Plots of various integrated FB asymmetries with the branching ratio of the inclusive 
mode (B XsT^t^) in the mSUGRA model. The mSUGRA parameters are: uiq = SOOGeV, 
^1/2 = 450GeV, A = 0. In these plots we have varied tan/3 in the range 25 < tan/3 < 50. 

analysis we will use a more generous bound [19]: 

2 X 10"^ < Br{B Xsj) < 4.5 x 10"^ 

which is in agreement with CLEO [20], ALEPH [21], BELLE [23] and BaBar [24] results. 

We have checked that the qualitative natiirc of our results doesn't change even if we use a 
much stronger bound on the branching ratio of i? ^ Xg^- As has already been mentioned 
in many earlier works [11] it is the purely dileptonic decay, Bg which is most 

affected by the introduction of the new set of operators. Therefore this decay would be 
one of the most promising modes to search for signatures of these new operators. The SM 
prediction of this mode is ~ 10~^ , which would be enhanced by many orders by these new 
operators. This process {Bg —>■ has not yet been observed in B-factories but CDF 

has given an upper limit on this branching ratio [25]. For our numerical analysis we have 
considered only that region of the mSUGRA parameter space which complies with the new 
CDF Run II bound Br{Bs -> ^+^-) < 9.5 x 10"^ [25] 

^^we wish to thank the referee for updating us on the new CDF result, in our first version of the manuscript 
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In Figure ^ we have shown the plots of the branching ratio of the process concerned 
as a function of tan/3 for various values of mg- In Figure ^ we have shown the variation of 
the Ann with the invariant mass of the dileptons (s) for various values of tanf3. As can be 
seen from the expression of ^jVA^ given in equation ( |3.1S| ) (and also given in Table |l|) the 
asymmetry vanishes within the SM, as within the SM the scalar and pseudo-scalar operators 
are negligibly small (~ m^mfe/m^). Therefore any observation of this asymmetry can be 
considered as a signal of new physics. In Figure ^ we have shown the values of various 
integrated FB asymmetries as a function of the branching ratio. In Figure |^ we have shown 
the contour plots Ann in the [mQ,mii2] plane for two different values of tanfi. 

As has been mentioned in many earlier works [2, 10-12, 14], the effect of the scalar (and 
pseudo-scalar) operators crucially depends on tan[3 and the mass of the Higgs. One could 
have used some relaxed kind of SUGRA framework also. By relaxed we mean relaxing 
the universality of the masses and coupling constants. Here we can relax the conditions 
of universality of scalar and/or gaugino masses. In the literature these models have been 
termed as rSUGRA (or relaxed SUGRA models) [18]. In these situations one can have 
sufficiently low allowed^^ values of the Higgs mass for large tanj3. In such situations the 
numerical values of the Wilson coefficients accompanying the scalar and pseudo-scalar op- 
erators could be high, and hence their effects on observables could be more profound. 
Similarly one could have used CMSSM models, where even the condition of correct elec- 
troweak symmetry breaking is relaxed. In these models one can have a very large variation 
in the Wilson coefficients Cq^ and Cq^- However, in this work we have confined ourselves 
to the mSUGRA model. 
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Figure 4: Contour plots of various values of averaged < Ann > (as given in plots) in the [mo, mi/2] 
plane for the mSUGRA model with A = and (a) tanP = 40, (b) tan/S = 45. 

So far no experimental search for the inclusive decay mode B XgT^r^ has been 
carried out. Experimentally this decay can be searched for by observing the /x produced 
from the r via the r — > fivi? process, but the muons produced in this cascade decay would 
be soft ones (so to be observed one has to observe muon pairs at low invariant mass). But 
there was an upper bound put on the branching ratio of this mode {B X^r+r") by 

we considered the older CDF bound [26] Br{Bs ^^'^^^~) < 2.0 x lO"*^ 

^^by allowed we mean that values are consistent with experimental constraints. 
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Grossman et al. [28] which was: 

Br{B XsT+T^) < 5%. 

As the r detection efficiency is low, such that at the SM level this decay mode might 
remain out of reach of the present generation B-factories (like BaBar, BELLE etc.), it has 
been well emphasized in many earlier works [3, 14, 28] that this decay can be enhanced by 
a couple of orders in magnitude in the presence of new physics and hence could be very 
useful in putting bounds on new physics models. 

To summarize the major conclusions of our numerical analysis: 

1. Almost all the polarized asymmetries show a large variation from their respective SM 
values over a large parameter space at large tan(5. 

2. The SM relationship that the magnitude of All is the same as A is violated in the 
presence of new operators. 

3. The observation of a non- vanishing value for Ann could be treated clearly as a signal 
of new physics. 

These polarized FB asymmetries also provide us with much needed observables, which 
could aid in the testing of the effective Hamiltonian's structure; especially the scalar and 
pseudo-scalar parts as can be see from the expression for Ann which is proportional to 
Cq^ and hence vanishes within the SM. 
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A. Input parameters 



TTT-T- 

a 



rrib = 4.8 GeV , = L4 GeV 
1.77 GeV , = 80.4 GeV , 

91.19 GeV , VtbVtl = 0.0385 , 
Gf = 1.17 X 10^5 GeV-2 
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